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Estimators for weak lensing observables such as shear and convergence generally have non-linear 
corrections, which, in principle, make weak lensing power spectra sensitive to primordial non- 
Gaussianity. In this paper, we quantitatively evaluate these contributions for weak lensing auto- and 
cross-correlation power spectra, and show that they are strongly suppressed by projection effects. 
This is a consequence of the central limit theorem, which suppresses departures from Gaussianity 
when the projection reaches over several correlation lengths of the density field, Lp ~ 55Mpc//i. 
Furthermore, the typical scales that contribute to projected bispectra are generally smaller than 
those that contribute to projected power spectra. Both of these effects are not specific to lensing, 
and thus affect the statistics of non-linear tracers (e.;;., peaks) of any projected density field. Thus, 
the clustering of biased tracers of the three-dimensional density field is generically more sensitive to 
non- Gaussianity than observables constructed from projected density fields. 
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I. INTRODUCTION 

In recent years, a renewed interest in the effects of 
primordial non-Gaussianity on the large-scale structure 
of the Universe has emerged, prompted by the significant 
effect on the bias of dark matter halos at large scales 
measured in N-body simulations [ll-Q- The observation 
of this effect in redshift surveys would be able to provide 
an independent confirmation of a possible detection of 
primordial non-Gaussianity from the anisotropies of the 
cosmic microwave background (CMB). Such a detection 
would open a completely new perspective on infiation 
and the high-energy physics of the early Universe [see, 
for instance, 0|- 

In fact, a specific type of non-Gaussian initial condi- 
tions, i.e. the local model of primordial non-Gaussianity 
d, [§], induces a strongly scale-dependent correction to 
the linear halo bias. This correction has been derived 
using several approaches, mostly based on the peak- 
background split formalism P, [2, llol-fTsI or on the statis- 
tics of high-peaks (also, see |l6t for a different per- 
spective). On the other hand, one can make the simple 
assumption of a local but nonlinear bias relation between 
the galaxy distribution and the underlying matter field 
■ When applied to the local model of non-Gaussianity, 
this yields the same scale-dependent correction as ob- 
tained in the former two approaches [l^, [l^ . 

The effect due to nonlinear local biasing shows an es- 
pecially close analogy with the case of non-linear weak 
lensing estimators we consider in this paper. A nonlin- 
ear but local galaxy bias relation can be expressed by the 
Taylor expansion [20| 



where 5g{x) and 5{x) represent, respectively, the galaxy 
and matter density contrasts and bi and 62 are the linear 
and quadratic bias parameters, which are assumed to be 
constants for a given galaxy sample. This leads to the 
following perturbative expansion for the galaxy two-point 
correlation function: 

{6g{Xl)5g{X2))^bl{5{Xl)5{X2)) 

+ ^&l&2 [{5{xi)5^{x2)) + {xi O X2)] 

+ •••, (2) 

where the second term on the r.h.s. depends on the mat- 
ter 3-point correlation function or bispectrum. For Gaus- 
sian initial conditions, the matter bispectrum induced 
by gravitational instability leads to negligible corrections 
in the two-point correlation on large scales. For non- 
Gaussian initial conditions of the local type, on the other 
hand, the correction becomes relevant on large scales. 

Since such behavior arises simply from the nonlinear 
relation between the observed distribution 5g and the 
matter distribution 5, we expect similar effects to show 
up for any other large-scale structure observables where 
analogous nonlincarities are present. 

Weak lensing measures the three components k, 71, 72 
of the tidal tensor of the lensing potential '0, 
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5g{x) = bi 5{x) + -b2 5\x) + 
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where derivatives are with respect to angular coordinates 
on the sky (we assume small angles, zero curvature, and a 
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flat sky limit throughout), x denotes comoving distance, 
while Xs is the distance to the source galaxy being lensed, 
and $_ = ('5 — $)/2is the lensing potential in conformal 
Newtonian, or longitudinal, gauge. 

While the two components of the shear 7,; can be mea- 
sured using galaxy ellipticities, the convergence k can be 
measured either through its effects on the number den- 
sity of galaxies (magnification bias), or on galaxy sizes 
and fluxes. In all of these cases, the estimators are not 
purely linear. This is not surprising since Eq. ([3]) is only 
a lowest order approximation to the lensing effect. Given 
that the shear 7 is a (two-dimensional) spin-2 field while 
K is a scalar, we can generally write 



7°'^^(0)=7(0) 



^obs 



c^k(6/)7(6/) + ■ 

{e)^K{e) + c^^K'{e) 



(5) 
(6) 



where the measured convergence and shear, 

and the dots indicate third- and higher order terms. For 
simplicity, we have assumed that all additive and mul- 
tiplicative biases have been removed so that to lowest 
order, k°^^ — k and 7°'^'' = 7. 

In the case of shear measurements from galaxy surveys, 
corrections to the leading order come from two sources: 
the fact that galaxy shapes estimate the reduced shear 
(7 = 7/(1 — k) rather than the shear 7 [2l| - |23| : and lens- 
ing bias [l^l, the fact that we preferentially select lensed 
galaxies in a flux- and size-limited source galaxy sample. 
Following the estimate of [2^, these two effects add up, 
yielding Cj roughly between 2 and 3. 

Examples of measurements of the convergence include 
using the number density of backgrou nd g alaxies via the 
magnification and size bias effects [H, [27f . and using the 
sizes and other measured characteristics of galaxies [28l - 
Isoj . All these estimators in fact measure the magnifica- 
tion /X, 
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(l-^)2-|7|2 



1 + 2k + 3k2 + |7|2 + ... (7) 



Hence, in the ideal case we expect c^i = 3/2 and Ck2 = 
1/2. These values will likely be modified in practice due 
to galaxy selection effects similar to those present in the 
shear. Note that it is often possible to vary the non- 
linearity coefficients , Ck experimentally, for example 
by applying different cuts on the source galaxy sample. 

In analogy with Eq. the two-point correlation of 
a non-linear estimator such as Eqs. (O-® receives cor- 
rections from three- and higher point functions of the 
underlying density field. This is easy to see e.g. for the 
shear two-point function (neglecting the spinor indices 
for simplicity): 



(7°b^(0)7°i^^(0')) = {l{0h{0')) 



(8) 



The leading correction term is a three-point function of 
shear and convergence, evaluated in the limit where two 



of the three vertices coincide. In [2J, |3l|, |321 , this con- 
tribution was investigated for non-Gaussianities from the 
non-linear evolution of the matter density. However, pri- 
mordial non-Gaussianity also modifies the shear and con- 
vergence power spectra in the same way. Throughout 
this paper, we shall focus on the local model of non- 
Gaussianity, where the bispectrum of primordial curva- 
ture perturbations, evaluated at last scattering, is given 
by 



B^{ki, fc2, fca) = 2fNi^[P^{ki)P^{k2) + perm. 



(9) 



Our conclusions however are broadly valid for any pri- 
mordial bispectrum shape. 

It will be useful to compare the impact of non- 
Gaussianity on weak lensing correlations with the analo- 
gous effects on the clustering of large-scale structure trac- 
ers such as galaxies or halos. As shown in [l[ , tracers with 
a linear (Eulerian) bias 61 acquire a scale-dependent bias 
correction in the presence of local non-Gaussianity given 

by 



A&i(fc; z) 



2hLibi - 1)4 
M{k,z) 



(10) 



where Sc = 1.686 is the critical density for the spheri- 
cal collapse in the flat matter dominated universe. The 
function A4{k,z), which relates the matter density fluc- 
tuations to the initial curvature perturbations as 5(k) = 
A^(fc, z)(/)(k), is given by 



M{k,z)^'^^D{z). 



(11) 



where T(fc) is the matter transfer function, Cp = 
3/2 fi^iJQ, and D{z) is the linear growth function nor- 
malized to the scale factor at last scattering (at which 
in our convention /nl is deflned), that is D{zis) = 

1/{1 + Zls). 

Throughout this paper, we shall also assume the Lim- 
ber approximation and work in the small angle limit, 
where the spherical harmonic transform is reduced to 
the two-dimensional Fourier transform. This very use- 
ful approximation does not significantly infiuence our re- 
sults [33 . We also assume that sources reside at a fixed 
redshift. Our fiducial cosmology is defined through the 
maximum likelihood cosmological parameters in Table 1 
("WMAP5+BA0+SN") of Komatsu et al. [p. 

We begin in Sec. |ll] by studying the impact of non- 
Gaussianity on shear and convergence power spectra. We 
then investigate the cross-correlation of shear and conver- 
gence with large-scale structure (LSS) tracers in Sec. IIIII 
Sec. [IV] extends the arguments to the general case of the 
clustering of peaks (or more generally nonlinear trac- 
ers) identified in projected density fields. We conclude 
in Sec. IV] 



3 



II. SHEAR AND CONVERGENCE POWER 
SPECTRA 

We first consider the impact of primordial non- 
Gaussianity on large-scale shear and convergence power 



spectra. By Fourier transforming Eq. ((8)) and the anal- 
ogous equation for the convergence, we obtain the ob- 
served shear and convergence power spectra: 



(27r) 
(2^)2 



[c«i + 2cK2Cos2(0f, - (f>i-i,)] B^iii, \e - £i\,£) 



(12) 
(13) 



Here, is the convergence bispectrum, and we have 
aligned so that = 0. Further, we have used that in 
Fourier space ^{£) = ^i{£) + i"f2{f-) = exp{2i(j)g)K{f.), and 
hence the leading order shear power spectrum is equal 
to the leading order convergence power spectrum Ci^{£). 
The convergence power spectrum and bispectrum are re- 
lated to the corresponding matter correlators P^, B^ 
via 

C.(£) = C^r^^^%fp„(£/x;x),(14) 
Jo X Xa\x) 

Jo X X^a'^ix) 

xB„ — ,— ,— ;X , 15 
\X X X J 

where Wl{Xs,x) = x/Xs{Xs - x) denotes the lensing 
kernel. 

In Eulerian Perturbation Theory (PT), the leading or- 
der expression ( "tree- level" ) of the matter bispectrum, 
valid at large scales, is given by the sum of a primordial 
component, i?/, present for non-Gaussian initial condi- 
tions, and a contribution Bq due to second-order cor- 
rections to the matter density induced by gravitational 
instability. We have 

B^^Bi + Bg + ... , (16) 

where the dots represent higher-order corrections in PT. 
The primordial component Bi is directly related to the 
primordial curvature bispectrum B^ via 

B/(fci, fca, fca; z) = M{ki, z)M{k2, z)M{k3, z) 

X B^(fcl,fc2,fc3), (17) 

while the non-linear component Bq is given by 

Baiki ,k2,k3;z)^ 2F2 (ki , k2)PL (fci ; z)Pl (fcs; z) 

+ (2 perm.) , (18) 

where PL{k; z) = ^/l{k,zYP^{k) is the linear matter 
power spectrum, and 

F2(ki, k2) = 5/7 + (l/2)ki • £2(^1/^2 + k^/ki) 

+ (2/7)(ki -£2)'. (19) 



Note that at small scales, additional perturbative correc- 
tions become relevant and, in general, it is not possible 
to separate out the purely primordial components in the 
matter bispectrum [3^, [3al ■ We shall return to this issue 
below. 

From the non-Gaussianity viewpoint, Eqs. (|12p - P^ 
come as no surprise: neglecting the distinction between 
7 and K for the moment. Eqs. ([S])-® say that the ob- 
served shear and convergence are biased estimates of the 
true 7, K with linear bias hi = 1 and quadratic bias 
parameters 62 = c^,Ck. Hence, Eqs. (fT2 |) -([T3 |) are analo- 
gous to the expression for the non-Gaussian halo power 
spectrum [l^, with two differences: first, we arc dealing 
with the projected density field k; second, the relation 
between shear and convergence leads to cosine factors in 
the integral. Note that on large scales, the bispectrum P^ 
is evaluated in the squeezed limit, since typically £ <^ £1. 
This is again similar to the halo clustering case. 

Fig.[l]shows the relative amplitude of the correction to 
the shear and convergence power spectra from a numer- 
ical evaluation of Eqs. P^ - ((T5)) when using the primor- 
dial bispectrum of the local type. We also show the tree- 
level bispectrum from non-linear evolution which con- 
tributes at the percent-level to C^oba and C^oba . Clearly, 
the contributions from primordial non-Gaussianity are 
strongly suppressed: they are always below 10^"* for the 
shear, and only reach 10"'^ at the very largest scales for 
the convergence. This is in stark contrast to the results 
for the halo power spectrum [l|, [lo| where order unity 
corrections are observed for this type of non-Gaussianity 
on large scales. 

In order to understand where this suppression comes 
from, we make an order of magnitude estimate of 
Eqs. (fT2l) - (fT3)) . We approximate the convergence power 
and bispectrum as 

C,{£)r.ClxxlPL{l/x:x), (20) 

and 

B.{£,j2,h)-Clx%BA^-^/-^/-^-i\ . (21) 

\X X X J 
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FIG. 1: Relative correction to the shear and convergence 
power spectra from primordial NG of the local type (red 
solid/dashed), and the tree-level matter bispectrum from non- 
linear evolution (blue long-dashed/dotted). Here, we have as- 
sumed sources located at = 2 and non-linear parameters 
as indicated in the figure (see Sec. |l]). 



Here, x is an effective lens distance, and xp, xb are di- 
mensionless geometrical factors of order unity which we 
leave unspecified for the moment. On large scales, the 
bispectrum is evaluated in the squeezed limit, /ci 3> fc. 
In the squeezed limit, Eq. (|9]) asymptotes approximately 
to 

BKfci,fc,|ki-k|) =2/NLA^(fc)X(fci)A^(Iki-k|) 
xP^(fc) [P^(fci) + P^(|ki-k|)] 



4/NLA^(fc)P4fc)^L(fcl) 



(22) 
(23) 



neglecting the third permutation which is suppressed. In 
the second approximate equality, we also set |ki— k| « k\. 
Further, for the moment all power spectra arc assumed 
to be evaluated at z(x)- After some manipulation, we 
arrive at the following estimates: 



AC^obs(^) «0, 
AC7„ob.(^) «8/nl 



Cpy^XB \xp 



where 



2 



(24) 



(25) 



(26) 



The cos 2(/>f J factor in Eq. ((T^ leads to a complete cance- 
lation of the effect on C^obs in this approximation, while 



FIG. 2: Effective lensing weight functions for the shear power 
spectrum (red dashed), and the shear bispectrum (blue solid). 



the phase factor in the second term of Eq. ([T5)) becomes 
unity. Further, note that in reality there is a high— ^ 
cut-off in Eq. (pS)) due to the resolution or pixel size of 
the shear survey. As long as this cut-off is less than ^ 
arcmin-scale, however, the quantitative results do not de- 
pend sensitively on the resolution. 

Using Eq. ([SO)) . Eq. ([25|) can be further simplified to 
yield 



8/] 



NL 



e^T{e/x) XB 



Xb 

Xp 



(27) 



where = c„i -f- Ck2- Note the appearance of an 
just as a factor of appears in the halo biasing in 
the local model of NG [l[. In fact, it is instructive to 
compare Eq. (|27| to a similar estimate for the angular 
power spectrum of some biased tracer "/i". Using the 
Limber approximation, we have 



dx Fh{xf 
X X 



Phh ( -;x 



(28) 



where Fh (x) is the selection function, normalized to unity 
in comoving distance. Let us now assume that the tracer 
is localized in a narrow redshift slice around a comoving 
distance x- Given that Phh{k) = [bi + Ab{k)]'^ Pmik) 
where A6 is given by Eq. (fTO)) in the local model of NG, 
we can approximately write the leading correction to Chh 
as 



2/i 



NL- 



Cpx' 



6i 'D{z)e^T{e/x)'' 



(29) 



5 



10- 



10- 



u 




10- 



FIG. 3: Relative amplitude of the non-linear lensing correc- 
tion from primordial non-Gaussianity to the shear and con- 
vergence power spectra as function of the source redshift Zs, 
for a fixed multipole I — 100. The other parameters are at 
the same values as for Fig. [1] 



where z = z{x). Assuming that % is a cosmological dis- 
tance, all factors in Eq. are in fact of order unity. 
Comparing Eq. (pS)) with Eq. (P7| shows that the sup- 
pression of the impact of non-Gaussianity on weak lensing 
power spectra comes from two sources: first, it is sup- 
pressed by CT^ 10~^. As we will show in Sec. IIV| this 
suppression factor is due to the projection of the density 
field over many correlation lengths, and thus ultimately 
a consequence of the central limit theorem. Second, 
the projection kernel for the bispectrum strongly prefers 
small lens distances, so that xb xp ~ 0.5. To illus- 
trate this, Fig. [2] shows the redshift weighting function of 
the shear bispectrum and power spectrum in comparison. 
Clearly, most of the contribution to the shear bispectrum 
comes from low redshifts, x/Xs ^ 0.2. In other words, 
for a survey with = 2 {xs ^ 3700 Mpc/ft.), the shear 
bispectrum hardly probes scales above 700 Mpc//i. Since 
the effect (at least from local NG) peaks on large scales, 
this is a severe limitation. 

Finally, choosing some numbers, cr^ = 10^"^, xb — 0.2, 
xp = 0.5, and setting T = 1, we have 



10" 



/nl 
100 



100 



(30) 



Comparing with Fig. [1] we see that the order-of- 
magnitude estimate predicts the right amplitude to 
within a factor of a few. We also see the behav- 
ior for AC^oba on large scales, and that AC^obs, is indeed 
strongly suppressed for small £ due to the cosine factor 



in Eq. (|12p . The restriction to small scales due to the 
lensing projection can be somewhat mitigated by going 
to larger source redshifts. Fig. |3] shows the evolution of 
the non-linear corrections with source redshift. While 
the corrections, especially AC^obs, increase with Zg, the 



corrections remain much smaller than 10 at 



100. 



A. Impact of non-linearities on small scales 

So far, we have only considered the leading order (tree- 
level) contribution to the matter bispectrum from pri- 
mordial non-Gaussianity. A simple estimate of the effect 
of non-linearities can be obtained directly from Eq. (|27p , 
by replacing cr^ with the non-linear variance of the con- 

5(7^ (estimated using the non- 



vergence. 



Since cr^ 



37|), 



linear matter power spectrum from halof it 
linearities are expected to increase AC°^''(€) by the same 
factor. We can obtain a somewhat more sophisticated es- 
timate by using the fact that, in the local model of non- 
Gaussianity, a long wavelength (linear) perturbation (pi 
acts to increase the variance of small-wavelength matter 
perturbations 5l,s [ll. IToj: 



{5U 



(dl,) (1 + 4/NL0O, 



(31) 



where all perturbations here are evaluated at early times, 
i.e. in Lagrangian space (indicated by subscript "L"). 
This is clearly reflected in the approximate squeezed-limit 
expression for Bm, Eq. (^5)) (since M{k)P^{k) = (Sicpi)). 

In other words, a given region in such a non-Gaussian 
Universe with a fixed value of /nl^; > is statistically 
equivalent to a region in a Gaussian Universe with a 
slightly higher amplitude of the primordial power spec- 
trum. Thus, we can model the change in the statistics of 
the late-time, non-linear small-scale modes (in Eulerian 
space) in the presence of a fixed large-scale mode as 



{SEAk)SEA-k)) 

dPnAk) 



X 1 + 4/NL0r 



din As 



(32) 



where PNL(fc) is the non- linear matter power spectrum 
in a ACDM cosmology with Gaussian initial condi- 
tions. Specifically, we estimate the contribution to the 
non-linear matter bispectrum due to primordial non- 
Gaussianity in the squeezed limit (ki ^ k) as 

BjMki,k,k') ^ 2h^M{k)M{ki)M{k') 



X P0(fc) 



M-^{k,) 



dPMk) 
a In A 



+M-'ik') 



,MAk') 



din As 



(33) 



with k' = ki — k. Since Pl As, this equation recov- 
ers Eq. when all fc-vectors are in the linear regime. 
We use the halof it prescription to evaluate the deriva- 
tives of Pnl- Fig. 0] shows the correction to the shear and 
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FIG. 4: Same as Fig. [TJ but with non-linear correc- 
tions included. For the corrections from primordial non- 
Gaussianity (red solid/dashed), we use the matter bispectrum 
from Eq. ((55)) with Pnl(/c) given by halof it [sj- For the cor- 
rections from non- linear evolution (blue long-dashed/dotted), 
we use the fitting formula from [Sg]. All other parameters as 
in Fig. m 



quantities given by 

= [6i(zl) + A5 (fc = ^/xl; zl)] 

X (1 + ^^)E£%lA^P„^ (^11^^- , (34) 

Here, we have again used the Limber approximation. Chk 
can be estimated either through convergence estimators 
or through shear, by using the relation between k and 7 
in Fourier space. 

The derivation of the nonlinear lensing corrections pro- 
ceeds in close analogy to Sec. [TTl and we obtain: 



■ cos 2<j)i^ 



■'^ ' (2^)2 



(35) 



Ck.2 



(2^)2 



cos2{6t. 



(36) 



Here, the halo-K-K bispectrum is defined through 

(/i(£i)k(^2)ac(£3)) = (27r)2fo(^i23)B,.««(^l,^2,^3). (37) 

In the Limber approximation (and assuming linear bias- 
ing), it is given in terms of the matter bispectrum by: 



convergence power spectra when using the matter bispec- 
trum from Eq. Note that for £ > 400 the squeezed 
limit is not a good assumption anymore and Eq. (|33p 
loses its validity. We see that AC^oba is indeed boosted by 
matter non-linearities by a factor of 8 on large scales, 
comparable to the simpler estimate using Eq. (PS)) . On 
the other hand, AC^obs is still suppressed on large scales 
and remains insignificant. Fig. 2] also shows the correc- 
tion purely from non-linear evolution for Gaussian initial 
conditions, calculated using the bispectrum fitting for- 
mula from [3^ combined with halof it (this is the same 
prescription as adopted in [H, HH [s^l ) ■ This contribution 
is clearly boosted as well by including non-linearities on 
small scales, so that it still dominates over the correction 
from primordial non-Gaussianity up to very large scales. 



III. GALAXY-GALAXY LENSING 

We now consider the cross-correlation between a large 
scale structure tracer h (such as galaxies or galaxy clus- 
ters) and weak lensing shear 7 and convergence n. Such 
cross correlations, called galaxy-galaxy lensing, probe the 
tracer- mass cross power spectrum. As derived in [33l [39| , 
the two-point angular power spectrum between tracer 
and lensing convergence is to leading order in the lensing 



Bh..{£l.£2,£3;ZL) - [CpWLiXs,XL)il + zl)Y — 

XL 



X [bi + Abi£i/xL;zL)] 
X B 



h il ^3 \ 

— , — , — -.zl]- 38 

XL XL XL 



Again, we assume sufficiently large scales so that i?,„ = 
Bi -I- Bq. Note that we have two contributions of order 
/nl, cx h\Bi and oc A6i_Bg, and a contribution of order 

Fig. [5] shows the numerical result for Eqs. (|55|) - (P51) . 
assuming lens galaxies at — 0.74 and source galaxies 
at zs = 2.0. The redshift of the source galaxies is chosen 
by requiring Xs = 2xl- We assume that the linear bias 
parameter of lens galaxies is 61 =2, and the amplitude 
of non-Gaussianity is given by /nl = 50. We adopt the 
same values for the non-linearity coefficients c^,c„ 1,0^2 
as in the last section. We again see that the correction to 
Ch'-i is suppressed, while C^k receives a correction that 
strongly increases towards large scales; in fact, at ^ < 10, 
the term oc A6 Bi (x £~^ becomes larger than the con- 
tribution from the tree level matter bispectrum for these 
parameter values. Following similar steps to Sec. [Hi we 
can obtain an order of magnitude estimate of the nonlin- 
ear correction to Chu from the local primordial bispec- 
trum Bi (for Ch-^ we again have a cosine factor which 
leads to a cancelation in the large-scale limit). Assuming 
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FIG. 5: Relative correction to the halo-shear and halo- 
convergence cross power spectra from primordial NG of the 
local type (red solid/dashed). The blue long-dashed/dotted 
lines show the correction from the tree-level matter bispec- 
trum from non-linear evolution in the absence of primordial 
non-Gaussianity. Here, we have assumed a lens redshift of 
ZL — 0.74 and sources at Zs = 2 so that Xs ~ 2xi- The 
dash-dotted line shows the correction to Ch~,,ChK. from the 
scale-dependent halo bias Eq. (|10p [via Eq. (|34|l ]. assuming 
hi = 2. 



the approximation of Eq. (|22p wc can write AC^j^obs as 



4 c« bi /] 



NL 



M{i/xL; zl) 



cl 



.XL 



xl 



(39) 



where is defined below Eq. (|27p . Hence the relative 
correction to the tracer-convergence cross-correlation be- 
comes 



ACfe^(^) 



CIxIWl{Xs.Xl)Lp 
PT{i/xL)D^{Q) 



(40) 



Note that the bias factor 6i -I- Ab{l/xL) drops out. In 
these expressions, a new length scale appears due to the 
projection. 



(2^ 



P™(fc,z = 0). 



(41) 



Lp can be seen as the one-dimensional coherence length of 
the density field at z = 0, in the sense that the variance of 



FIG. 6: Relative amplitude of the non-linear lensing correc- 
tion from primordial non-Gaussianity to the halo-shear and 
halo-convergence cross-correlation, as function of the lens red- 
shift Zl, for a fixed multipole I — 100. The other parameters 
are the same as for Fig. O 



the density field projected along a thin slab of thickness 
Ax at redshift z is given by 



2 

"5 proj 



D{z) 



1 2 



D{0) 



Lp 
Ax' 



(42) 



In our fiducial cosmology, Lp = 54.6 Mpc/h. 

Since on large scales T{£/xl) ~ 1, Eq. (^0]) again re- 
covers the scaling of the non-linear correction, simi- 
lar to what was found for the weak lensing power spec- 
trum. Compared to the change in Chni^) due to the 
scale-dependent halo bias in Eq. ((M)) (green dot-dashed 
line in Fig. which is exactly 1/2 the effect on the 
halo angular power spectrum Eq. (|29p . the effect of the 
non-linear lensing correction is suppressed by a factor of 



CpWl{xl,Xs)Lp 



LpHo 



0.02, 



(43) 



where we have assumed that is of order the horizon 
scale. While this is a significant reduction, the projection 
effect docs not suppress the correction in galaxy-galaxy 
lensing quite as severely as it does for the shear power 
spectrum. This is mainly because in the cross-correlation 
with biased tracers the dominant contribution arises at 
Zl, rather than at low redshifts as in the lensing auto- 
correlation. Comparing with Fig. [5j the order of magni- 
tude prediction again matches quite well. 

Fig.[n]shows the amplitude of the corrections to AC^k, 
AChj as a function of source redshift. Here, we have 
chosen the lens redshift such that Xs = 2xl in each case; 
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thus maximizing the lensing kernel Wl- The scahng in 
redshift is very similar to that seen in the weak lensing 
power spectra, Fig. [3l 



IV. GENERAL STATEMENTS ON PEAK 
CLUSTERING IN PROJECTED FIELDS 

The results derived in the last two sections allow us to 
make some interesting and general statements on the ef- 
fectiveness of using the clustering of peaks (or, more gen- 
erally, non-linear tracers) identified in two-dimensional, 
projected density fields. Applications of this include 
shear peaks as well as peaks identified in diffuse back- 
ground maps. 

Consider a general projected density field 



F\{x) is peaked at some distance x- We then obtain the 
following estimate for the impact of non-Gaussianity on 
the clustering of peaks in the projected field A: 



NL 



/NL 



1 



C2 Cp^ 



D{-z) 



ciPT{e/x)D^O) ' 'Axx 



(50) 



m 



dxFxix)SixO;x), 



(44) 



where F\{x) is a filter function normalized to unity in 
X- Using the Limber approximation, the angular power 
spectrum of A is then straightforwardly obtained as 



XX \x 



(45) 



Here, z = z{x) and Axa = [/ dxF^{x)]~^ is the effective 
width of the projection kernel. In the second line, we 
have introduced the correlation length Lp via Eq. (|42p . 

This result summarizes the different expressions for 
non-Gaussian corrections to angular power spectra de- 
rived in this paper: the first line of Eq. ([SO]) clearly 
shows a strong similarity to the expression for halo clus- 
tering, Eq. (P^. upon identifying C2(j| — Ci5c from 
Eq. Inserting the appropriate value for the width 

Axl = [Cp(l + zl)Wl]~^ of the lensing kernel (at fixed 
source redshift), setting ci — 1, C2 = c„, and dividing by 
a factor of 2 we also recover the expression for galaxy- 
galaxy lensing Eq. (|40p (the factor of 2 comes in since 
Eq. ([50)) is for the auto-correlation while Eq. ([M)) is for 
a cross-correlation). Finally, we see that the quantity 
(T^ Lp/AxL is suppressed by the same projection ef- 



Now assume we identify peaks in A; that is, we apply 
some nonlinear transformation so that the peak density 
is perturbativcly given by 



Apk(0) = ciA(0) + |A2(0) 



(46) 



The simplest example is thresholding, i.e. Apk(^) — 
Q{X{9) — Ac), where 8 is the Heavisidc function. If 
Ac 3> (Ta, f a being the r.m.s. fluctuation of A given by 



dxFl(x) f 



X X 



ci, C2 attain the well-known values [ij 



Ac 



feet. 

In summary, the clustering of peaks identified in the 
projected density field is suppressed relative to the angu- 
lar clustering of peaks identified in the three-dimensional 
density field [Eq. (HH)] by a factor of D'^{z) Lp/Axx- In 
fact, if the kernel is broad, then the contributions to the 
bispectrum B\ arc dominated by low redshifts, as we 
have seen in Sec. [ill which further suppresses ACx^^^ be- 
yond Eq. ([50]) . Note that the Limber approximation and 
hence Eq. (|50p break down if Ax a ^ Lp, i.e. for very 
narrow kernels. 

Thus, unless the line-of-sight projection Eq. (j^J) 
-^Tn(^/x;x)j (47) (jan somehow be restricted to a range Axa ^ 

D^{z)Lp « 54.6 D'^{z)Mpc/h, the impact of primordial 
non-Gaussianity on the clustering of peaks identified in 
a projected density field is strongly suppressed. 



2 ' 



C2 



(48) 



The following considerations are completely independent 
of the precise peak definition and value of C2, however. 
Analogous to our derivations in Sec. |ll] and IIIIl the an- 
gular power spectrum of peaks can be written as 

CxM)^clCx{e) + ACxM), 



ACxM)^c^c2 



Bx{t,iiAh-i\)- (49) 



The second term can in principle be used to probe pri- 
mordial non-Gaussianity. We can now use the Limber ap- 
proximation, and in the large scale limit make the same 
set of approximations as in the previous sections. As 
a further simplification, we will assume that the kernel 



V. CONCLUSION 

Observations of the large-scale structure in the Uni- 
verse offer very promising opportunities for probing the 
initial conditions of the Universe, complementing the 
searches for deviations from Gaussianity in the temper- 
ature anisotropics of the cosmic microwave background. 
Weak lensing is one of the most powerful probes of large- 
scale structure, as it directly probes the underlying mat- 
ter distribution, thus circumventing many of the uncer- 
tainties in the tracer-mass relation. 

In this paper, we have shown that shear and con- 
vergence power spectra are in principle subject to cor- 
rections from primordial non-Gaussianity, since weak 
lensing estimators generally have non-linear corrections. 
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However, projection effects severely reduce the impact 
of primordial non-Gaussianity, and the effects are gener- 
ally much smaller than those from nonlinear gravitational 
evolution. This also holds when approximately including 
the effect of non-linear growth on small scales. We also 
investigate the same effects for galaxy-galaxy lensing, i.e. 
the cross-correlation of shear and convergence with some 
large-scale structure tracer. In this case, the suppression 
is somewhat mitigated, but still significant. 

Finally, in Sec. IIVI we provide a general argument 
why this suppression is generic to searching for non- 
Gaussianity in two-dimensional projected density fields. 
If the projection occurs over a longer line-of-sight dis- 
tance than the one-dimensional coherence length of the 
density field, Lp k, 54.6 Mpc//i (at z = 0), the pro- 
jected density field is more Gaussian than the underlying 
3D density field, a consequence of the central limit theo- 
rem [ioj . Thus, the effect of primordial non-Gaussianity 
on the clustering of peaks identified in such a projected 
density field is suppressed by ~ Lp/Ax, where Ax is 
the width of the projection kernel. This applies to peaks 



identified in weak lensing shear maps (where Ax ~ Hq) 
as well as in diffuse backgrounds — unless a definite con- 
nection between such peaks and dark matter halos in 
the 3D density field can be made. Thus, we expect 
that the statistics of tracers of the 3D density field will 
generally provide tighter constraints on primordial non- 
Gaussianity than those of projected fields. 
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